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' Abstract 
> . 

— ^ We show that the phase transition previously observed in dynamical tri- 

' angulation models of quantum gravity can be understood as being due to the 

I creation of a singular link. The transition between singular and non-singular 

■ geometries as the gravitational coupling is varied appears to be first order. 

, Dynamical triangulations (DT) furnish a powerful approach to the problem 

of defining and studying a non-perturbative theory of quantum gravity. Simply put, 
^ ' the functional integral over (euclidean) four-geometries is defined as some scaling 

^ . limit of a sum over abstract simplicial manifolds. In two dimensions this approach 

D I has been very successful [|| . In addition to the calculation of gravitationally dressed 

' anomalous dimensions, the DT approach, being discrete, allows for the use of non- 

perturbative methods such as computer simulation. In two dimensions this has 
^ I yielded new results; the measurement of fractal dimensions characterizing the quan- 

^ tum geometry |^, power law behavior of matter field correlators on geodesic paths 

1^] and insight into the problem of formulating a renormalization group for quantum 
gravity [|]. 

The observation of a two phase structure for the four dimensional models was 
noticed early on |^ . This led to the speculation that a continuum theory of gravity 
could be constructed by taking an appropriate scaling limit near the phase transition. 
However, in |^], evidence was presented that the transition was discontinuous. Such 
a scenario would rule out arguably the simplest possibility - that a continuum theory 
could be obtained by tuning the theory to some critical (gravitational) coupling. 

Independent of these issues it was noted in that the crumpled phase of this 
model possessed a singular structure - that the sum over triangulations was domi- 
nated by those with a single link common to a very large number of four-simplices. 
The endpoints of this link, termed singular vertices, were shared by a number of 
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four-simplices which diverged hnearly with the total triangulation volumeQ. 

In this note we have studied how this singular structure changes as the grav- 
itational coupling is varied. We find evidence that the transition is driven by the 
creation of singular vertices as the coupling is lowered, singular vertices being absent 
in the branched polymer (large coupling) phase. We support our observations by 
appealing to a mean-field argument due to Bialas at al Q . 

The model we have investigated is defined by the grand canonical partition 
function 

Zgc = e-''4^-*+^«^° (1) 

T 

The coupling kq plays the role of an (inverse) bare gravitational constant 
conjugate to a discrete analog of the integrated curvature - the number of vertices 
Nq. In practice we use K4 (the bare cosmological constant) to fix the volume (number 
of four-simplices A'4) and consider a canonical partition function 

Zc(Ko,V^) = Ee""'^° (2) 

Tv 

The class of triangulations Ty with volume V is further restricted to be those 
with the global topology of the four-sphere together with a local manifold constraint 
- that the vicinity of any point should be homeomorphic to a four-ball. 

We use a Monte Carlo algorithm, described in to sample the dominant 
contributions to Zq- Fig. |l] shows a plot of the local volume of the most singular 
vertex < > (i.e the number of simplices sharing that vertex) against coupling kq 
for lattices with volume 8K, 16K and 32K. 

Deep in the crumpled phase it is clear that < ujq > scales linearly with 
system size as expected. For sufficiently large coupling kq it is small, the ratio 
< uq > /V approaching zero as y — > 00. The behavior of loq for intermediate 
coupling is somewhat complex - on the 8K lattice it falls rather slowly for kq < 1, 
kinks for kq ^ 1, falls more rapidly until ^ 1.8 and then plateaus until ~ 2. 
We will call the regime from kq ~ 1 to ~ 2 the mixed region. For couplings larger 
than this < ujq > rapidly approaches its asymptotic value. A similar picture is seen 
at 16K - the volume ujq falls slowly up to kq ~ 2 after which it falls quickly to a 
small plateau for kq ~ 2.2. For couplings larger than this the value of ujq quickly 
approaches its limiting value. Thus the behavior is similar to the 8K data except 
that the mixed region is now much smaller. The arrow on the plot indicates the 
approximate position of the usual phase transition as revealed by measurements of 
the vertex susceptibility y (^{Nq) — [Nq)"^^ for large volumes. It appears that the 
rightmost boundary of the mixed region coincides with the phase transition. At 
32K the mixed region has narrowed to the point where it can no longer be resolved 

^Singular vertices were first noticed in p| 
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Figure 1: Singular vertex volume vs kq 

and the data is best interpolated with a curve that undergoes a rapid variation in 
its gradient close to the phase transition. 

Further insight can be gained by looking at the fluctuations in this local 
volume Xo = y (("^o) ~ ('^o)^)- In Fig- @ it is clear that the data for small volume 
lies on a curve with two peaks - corresponding to the boundaries of the mixed 
region. These peaks appear to merge on the phase transition point as the volume 
is increased. For 32K these two peaks have already coalesced into one. 

The explanation for these finite size effects can be understood by plotting 
the maximal link local volume < ui > (fig. |3|). The left edge of the mixed region 
KQ(y) matches an abrupt change in < ui > as kq is increased - indeed by looking 
in detail at the triangulations we find for couplings kq in the mixed region the 
link between the singular vertices has broken - there is no identifiable singular 
link - merely a gas of a variable number of remnant singular vertices (one, two 
or more vertices with large local volumes). The rightmost boundary of the mixed 
region Kq{V) then corresponds to the final disappearance of the remaining singular 
vertices. This conclusion is strengthened by looking at the conjugate susceptibility 
Xi = y (("^i) ~ ("^i)^) (fig- @)- The peak associated with the breakup of the 
singular link appears to merge with the peak associated with the disappearance of 
singular vertices as the volume gets larger. 
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Figure 2: Fluctuations in singular vertex volume xo 



To summarize, our study of the singular structure indicates that at finite 
volume there are two pseudo-critical couplings Kq and Kq associated to the creation 
of singular links and vertices respectively. These appear to lock into a single critical 
point Kq as y — s- oc0. Furthermore, the gradient of the < ujq > curves in the tran- 
sition region increases with volume V. This is consistent with a finite discontinuity 
oi < lvq > in the infinite volume limit. Indeed if we interpret < lvq > /V as an 
order parameter then we can read off the order of the transition from the scaling of 
its corresponding susceptibility xo- The numerical data is certainly consistent with 
a linear scaling of the peak in xo with volume, characteristic of a first order phase 
transition. 

If we run simulations close to Kq for large volumes we see signs of meta- 
stability ~ the system tunnels back and forth between two distinct states - one is 
branched polymer-like and contains no singular vertex and the other contains one 
or more remnant singular vertices. This type of behavior has been seen before and 
is a strong indicator of a discontinuous phase transition. 

The correlation between the action Nq and singular vertex volume ujq can 
be seen easily in fig. ^ where time series of the two quantities are plotted. There 

^Similar results were obtained independently by B. Petersson and J. Tabaczek 
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Figure 3: Singular link volume vs kq 

is clearly a perfect (anti)correlation - fluctuations to larger vertex numbers are ac- 
companied by precisely matching decreases in the singular vertex volume. This is 
yet another indicator that the physics of the phase transition corresponds exactly 
to the fluctuations in the singular structure. 

It is possible to gain further insight into the nature of this transition by 
making a couple of assumptions. Consider the microcanonical partition function 
Zmc (^i-^o) gotten by fixing the number of vertices. Following let us define a 
local entropy Si associated with a vertex i. This is just the number of ways of gluing 
together the simplices in its local volume element - hence it is a function of that 
local volume Sj = Si{LOo{i)). We will first assume that Zmc can be approximated by 
enumerating all possible ways in which these local entropies can be assigned^. This 
is a drastic assumption but we are encouraged by the results of our simulations which 
support the notion that the transition is driven by fluctuations in the singular vertex 
structure. Secondly we will assume that for sufficiently small vertex density Nq/V 
a mean field approximation can be used in which the local entropies are treated as 

^ In the light of the results derived in [11], we should probably regard this as an effective theory 
vahd for understanding some of the properties of the phase transition - it seems that the leading 
behavior of the partition function can be derived from a partitioning of simphces amongst the 
(non-singular) {d — 2)-simplices - the effect of the singular structure being sub-leading. 
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Figure 4: Fluctuations in singular link volume xi 



independent. Thus we are led to an ansatz for the partition function of the form 



ZMciV,No)= l[siiu;iii))5iY.iOo{ 



cV 



(3) 



The final delta function reflects the fact that at least one global constraint 
restricts the local volumes uJo{i) ~ that the sum of these volumes is just proportional 
to the volume V. (c the constant of proportionality is equal to five). As argued 
in for Si {u)o(i)) we should take the number of triangulations of the three sphere 
bounding the vertex i. The precise functional form of Sj is unknown but previous 



numerical simulations |12] and an analytic calculation of the exponential bound IT 
lead to the form 



•Si {n) ~ exp {—an'^)e 



bn 



(4) 



The constant a ~ 0.6 — 7. Using the results of Q one can write the partition 
function of this model in the thermodynamic limit as 



Zmc = ^ exp ((- log X + b)cV + InF (A)iVo) 
zvr 



(5) 
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Figure 5: Correlation of vertex number Nq with singular vertex volume 

where the parameter A is the solution of the equation 

XF' (A) 

and 



(6) 



F{X) = Y: s' {u) a- (7) 

UJ 

The prime indicates that the leading exponential piece in st has been trivially 
factored out (it can be seen in the leading term in eqn. |5|) and so F depends only 
on the sub-exponential pieces in Sj. 

The model exhibits a two phase structure; for densities p = cV/Nq with 

F' (1) 

p > Pc = the parameter A = 1 and the system is in a collapsed phase where 
vertices of small order behave independently and the global constraint is satisfied by 
a small number of vertices which are shared by a number of simplexes on the order 
of the volume. This regime is identified with the crumpled phase of the dynamical 
triangulation model. Conversely for p < pc the parameter A varies between zero 
and one and the free energy varies continuously with p. The distribution of vertex 
orders then behaves as 
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p(n) ~L^(n)e-"'°s^ (8) 

This phase is then identified with the branched polymer phase in the dynam- 
ical triangulation model. As p (and hence A'o) is varied (by varying the coupling kq) 
the system moves between these two phases. This phase structure is seen for any 
value of the parameter a < 1. Indeed any sub-exponential or power law behavior 
will suffice. 

Our simulations are run in the canonical ensemble where a sum over vertices 
is taken 

Zc = ^ Zmc exp KoiVo (9) 
If we use the mean-field approximation for Zmc this integration can be 



done 1 14 1 . The model again exhibits a two phase structure with a first order phase 
transition separating the crumpled phase from the branched polymer phase. This 
then is in agreement with the results of our simulations. 

In our mean field arguments we have only considered the vertex sector of the 
model - clearly the links can be considered as possessing a local entropy too. In 
this case this entropy would be equal to the number of triangulations of the two- 
sphere dual to the link. Similar arguments would then lead to the prediction of a 
two phase structure for the links - a collapsed phase with singular links and a fiuid 
phase in which the link distribution behaves exponentially with link local volume. 
The coupling kq at which this transition occurred would seemingly be independent 
of the coupling at which singular vertices appear. Although our simulations are 
consistent with this scenario at small volume it appears that these two transition 
points merge in the infinite volume limit. Of course, the vertex and link sectors are 
not truly independent and it is presumably the presence of these residual constraints 
that is responsible for the coincidence of the two transitions. 

In summary, we have reported numerical results which strongly support the 
idea that the phase transition observed in 4d simplicial gravity is associated with 
the creation of singular geometries. These singular structures, composed of a single 
singular link between two singular vertices, dominate at strong coupling. As the 
coupling ^ is decreased they eventually disappear. This process appears to be 
discontinuous and is the origin of the first order nature of the transition. Simple 
mean field arguments support these conclusions. 

This work was supported in part by NSF Grant PHY-9503371 and DOE 
Grant DE-FG02-85ER40237. We acknowledge useful conversations with B. Peters- 
son. 
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